Microscopic mechanism for the 1/8 magnetization plateau in SrCu 2 (BOa)2 
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The frustrated quantum magnet SrCu2(B03)2 shows a remarkably rich phase diagram in an 
external magnetic field including a sequence of magnetization plateaux. The by far experimentally 
most studied and most prominent magnetization plateau is the 1/8 plateau. Theoretically, one 
expects that this material is well described by the Shastry- Sutherland model. But recent microscopic 
calculations indicate that the 1/8 plateau is energetically not favored. Here we report on a very 
simple microscopic mechanism which naturally leads to a 1/8 plateau for realistic values of the 
magnetic exchange constants. We show that the 1/8 plateau with a diamond unit cell benefits 
most compared to other plateau structures from quantum fluctuations which to a large part are 
induced by Dzyaloshinskii-Moriya interactions. Physically, such couplings result in kinetic terms in 
an effective hardcore boson description leading to a renormalization of the energy of the different 
plateaux structures which we treat in this work on the mean- field level. The stability of the resulting 
plateaux are discussed. Furthermore, our results indicate a series of stripe structures above 1/8 and 
a stable magnetization plateau at 1/6. Most qualitative aspects of our microscopic theory agree well 
with a recently formulated phenomenological theory for the experimental data of SrCu2(B03)2- 
Interestingly, our calculations point to a rather large ratio of the magnetic couplings in the Shastry- 
Sutherland model such that non-perturbative effects become essential for the understanding of the 
frustrated quantum magnet SrCu2(BOs)2. 

PACS numbers: 05.30.Jp, 03.75.Kk, 03.75.Lm, 03.75.Hh 



I. INTRODUCTION 

Strongly frustrated quantum magnets in an external 
field are fascinating systems because the interplay be- 
tween interactions and kinetics can lead to very rich 
phase diagrams. The magnetization can be described as 
a gas of bosonic particles whose density is controlled by 
the external magnetic field, and because frustration typ- 
ically reduces the kinetic energy, Mott insulating phases 
(corresponding to magnetization plateaux) 1-3 , superfluid 
or even supersolid phases have been predicted to occur 4 ' 5 . 
The experimental observation of these phases is an on- 
going challenge. A major player in the field is the lay- 
ered copper oxide SrCu2(B03)2, in which several magne- 
tization plateaux have been observed 6-11 . However, the 
definitive sequence of plateaux and the presence of su- 
persolid phases remain open issues that call for further 
experimental and theoretical investigation. 

The magnetization of SrCu 2 (B03)2 is expected to be 
described by the two-dimensional spin- 1/2 antiferromag- 
netic Heisenberg model known as the Shastry-Sutherland 
model 12 in a magnetic field 

H = J'Y f S i .S j + jY f S i - S, -B^S? , 

<ij> •Ci,j3> i 

with J'/ J ~ 0.65, where the <Ci,j^> bonds build an array 
of orthogonal dimers while the <i,j > bonds are best seen 
as inter-dimer couplings (see Fig. 1). For J'/ J smaller 
than the phase transition point at ~ 0.7 13 ~ 16 , the ground 
state of the model is exactly given by the product of 
dimer singlets, and the magnetization process can be de- 
scribed in terms of hardcore bosons which represent po- 
larized triplons ji 1 ) = | ft) on the dimers interacting and 



moving on an effective square lattice 4 ' 17 ' 18 . 

All theoretical approaches agree on the presence 
of magnetization plateaux at 1/3 and 1/2 4 ' 16 ' 17,19-21 , 
in agreement with experiments 6 ' 8,11 . Additionally, a 
plateau at 2/5 has been recently proposed in Ref. 11. 
However, the structure below 1/3 is rather controversial. 
On the experimental side, the original pulsed field data 
have only detected two anomalies interpreted as plateaux 
at 1/8 and 1/4 6 , but the presence of additional phase 
transitions and of a broken translations! symmetry above 
the 1/8 plateau has been established by recent torque 
and NMR measurements up to 31 T 9 ' 10 . The possibility 
of additional plateaux has been pointed out by Sebastian 
et al 8 , who have interpreted their high- field torque mea- 
surements as evidence for plateaux at 1/q with 2 < q < 9 
and at 2/9. 

On the theoretical side, the situation is not set- 
tled either. The finite clusters available to exact 
diagonalizations prevent reliable predictions for high- 
commcnsurability plateaux, and the accuracy of the 
Chern-Simons mean-field approach initiated by Misguich 
et al. 21 and recently used by Sebastian et al. s to ex- 
plain additional plateaux is hard to access. The essen- 
tial difficulty lies in the fact that, since plateaux come 
from repulsive interactions between triplons, an accurate 
determination of the low-density, high-commensurability 
plateaux requires a precise knowledge of the long-range 
part of the interaction. 

A quantitative step in calculating this long-range part 
of the interaction has been achieved by high-order series 
expansions 18 as well as by contractor renormalization 22 . 
In Ref. 18 an effective low-energy hardcore boson model 
for the Shastry-Sutherland model in a magnetic field is 
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FIG. 1: (Left) Illustration of the Shastry-Sutherland lattice 
and of the two-body interactions Vshjfij+s. Thick solid lines 
(dot-dashed lines) correspond to the magnetic exchange cou- 
pling J (J'). The two-body interaction labelled by Vs are 
defined as the density-density interaction between the thick 
dirtier labeled as j and the dimer labeled Vs. (Right) The 
hopping amplitudes t\ and ti are illustrated. 



derived. The classical solution of the effective model 
(which is expected to work very well at low densities 
n < 1/6 and not too large values of the perturbation 
J' I J R3 0.5) gives a sequence of magnetization plateaux 
at 1/9, 2/15, and 1/6. Interestingly, the experimentally 
most prominent and most studied 1/8 plateau is not fa- 
vored energetically 23 . 

Recently, a phenomenological theory based on inter- 
preting boron and copper NMR data of SrCu2(B03)2 
in high magnetic fields reveals an even more complex 
magnetization process 24 . Evidences for a sequence of 
plateaux at densities 1/8, 2/15, and 1/6 is found. Fur- 
thermore, the regime between 2/15 and 1/6 is interpreted 
as an infinite hierarchy of stripe structures. Interest- 
ingly, the structure of the different plateaux is in dis- 
agreement with the proposed plateau structures in the 
classical limit 18 . 

In this work, we are aiming at a microscopic theory 
which we want to compare with the phenomenological 
findings. To this end, we treat the quantum fluctuations 
of the pure Shastry-Sutherland model present at low den- 
sities and we extend the effective low-energy model de- 
rived by perturbative continuous unitary transformations 
(pCUTs) 18 by the dominant effect of additional magnetic 
couplings. Most importantly, we find that the inclusion 
of Dzyaloshinskii-Moriya (DM) interactions results in a 
natural microscopic mechanism for the stability of a 1/8 
plateau with a diamond unit cell. Physically, such un- 
frustrated terms give rise to kinetic processes in the ef- 
fective description and therefore they introduce quanum 
fluctuations in the classical plateau structures. Addi- 
tionally, it is shown that our theory predicts a stable 1/6 
plateau whose structure is only consistent with the phe- 
nomenological theory if the coupling ratio J'/ J is rather 
large. This is not a consequence of the induced quantum 



fluctuations. It results from the fact that the two-body 
density-density interactions do not respect anymore the 
perturbative hierarchy close to the phase transition point 
and therefore novel structures are stabilized. Finally, our 
theory predicts a series of stripe structures between 1 /8 
and 1/6. Altogether, our theory shares many similari- 
ties with the recent findings for SrCu 2 (603)2 in a mag- 
netic field. Nevertheless, several discrepancies between 
the microscopic and the phenomenological theory still 
remain. This is either due to the limited access to the 
most complicated non-perturbative regime of the stud- 
ied microscopic model or due to the fact that the physics 
of SrCu2 (1303)2 in a magnetic field depends on further 
subtle details of the material. 

The paper is organized as follows. In Sect. II, we ex- 
plain the origin of the classical plateaux in the effective 
model and we introduce the relevant plateau structures. 
In Sect. Ill, the additional couplings incorporated in this 
work are introduced in detail and their most important 
physical effect is illustrated. Additionally, the details of 
the mean- field approach are given Sect. IV. We finally 
present our results in Sect. V and we summarize the ma- 
jor findings in Sect. VI. 



II. CLASSICAL PLATEAUX 
A. pCUT+CA 

The pCUT transforms the Shastry-Sutherland model 
into an effective model conserving the number of 
triplons 18 . The relevant processes for the physics in a 
finite magnetic field have maximum total spin and total 
S z . This is true as long as bound states are not essential 
at low energies 25 which we assume in the following. The 
effective Hamiltonian obtained by the pCUT takes then 
the form of an interacting hardcore boson model where 
the amplitudes of the hardcore boson model are given as 
a high-order series expansion in J'/J. 

Typically, any kind of kinetic and interaction processes 
are present in such an effective model. The effective 
Hamiltonian H c s is by no means simpler than the orig- 
inal one in general, but it is in the limit of small den- 
sity. In Ref. 18 all terms with up to three creation 
and annihilation operators and all four-body interactions 
that first appear up to order < 8 have been kept. But 
in the small density limit n < 1/6, the magnetization 
is not affected by the three-particle and four-particle 
interactions 18 . Furthermore, in that limit the kinetic 
terms are very small, and they can be considered as a per- 
turbation of the interaction part. It is thus appropriate 
to use a Hartree approximation in which the variational 
ground state is a product of local boson wave-functions 
since this approximation becomes exact in the limit of 
vanishing kinetic energy. 

This Hartree approximation is most simply imple- 
mented by mapping the effective model onto a spin 1/2 
model using the Matsubara-Matsuda representation 18 ' 26 
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of hardcore bosons. In the spin language it then trans- 
lates into the classical approximation (CA) where the 
spins are treated as classical vectors of length 1/2. 

It has been found in Ref. 18 that only magnetiza- 
tion plateaux are realized in the classical limit (except 
a tiny supcrfiuid in the dilute limit) although superfluid 
or supersolid phases can be described within the CA. 
The magnetization plateaux correspond to Wigner crys- 
tals where triplons are frozen in a periodic fashion in the 
ground state breaking the discrete translational symme- 
try of the problem. The classical energy of all plateaux 
at low densities n < 1/6 can be quantitatively captured 
by the following effective Hamiltonian 

3 $ 3 

where the sum over j runs over all the sites of the effective 
square lattice formed by the dimers. The first term repre- 
sents the chemical potential (fi — (J>o) °f hard-core bosons 
where /i corresponds to the external magnetic field B and 
where fio originates from the Shastry-Sutherland model. 
We have calculated (iq up to order 17 in J' / J (see ap- 
pendix VIII). The second term denotes the two-triplon 
density-density interactions Vg. Note that the sum over 
6 contains all interaction terms illustrated in Fig. 1 plus 
their three symmetric counterparts. The classical energy 
of a given magnetization plateau depends therefore only 
on the arrangement of the local densities n = b'b . 

All two-body interactions V$ which first appear at or- 
der less or equal to 10 have been calculated up to order 
15 (except V\ which has been determined up to order 
14). The evolution with J'/J of the two-body interac- 
tions defined in Fig. 1 is depicted in Fig. 2. At small 
J'/ J, interactions beyond V4 are small and may be ne- 
glected, but for larger J'/J the higher order terms V3, 
V5, and V-j (appearing at order 6) become important and 
contribute to the formation of low-density plateaux. For 
these terms, the bare series and the dlogPade extrapo- 
lations are basically indistinguishable below J'/J = 0.6. 
Beyond that value, various extrapolations still give con- 
sistent results for the two-body interactions. Interest- 
ingly, we find that the perturbative hierarchy becomes 
invalid for J'/J > 0.65, i.e. one observes V3 ~ V\ and 
V3 ~ V4 although these interactions do not originate 
from the same perturbative order. In the following we 
deduce several indications that the frustrated quantum 
magnet SrCri2 (1303)2 is most likely situated in this most 
challenging but also most interesting non-perturbative 
regime. 

B. Plateaux structures and energies 

Representative results for the classical limit at J'/J = 
0.5 and J'/J = 0.68 are shown in Fig. 3. Physically, the 
realization of all relevant plateaux in the classical limit 
at low magnetizations follow the same guiding principle 











































/ %X 








/ft 


0.04 

// 8 / 


















/ 8 


5 0.7 



0.1 ~02 03 OA 0^5 0.6 0.7 



J'/J 



FIG. 2: Amplitude Vs of the extrapolated two-body interac- 
tions as a function of J'/J. Inset: Different dlogPade extrap- 
olants (solid lines) as well as the bare series (dashed lines) for 
V3 and V5. The displayed curves correspond to the data given 
in Ref. 18. 



for J'/J < 0.67. Triplons are placed such that only two- 
body density-density interactions Vg which start in order 
6 or higher have to be paid. Interactions starting in order 
6 are V 3 ', V 5 , and V 7 (see Figs. 1 and 2). Note that 
only the classical plateau at 1/9 does not pay any order- 
6 interaction. The explicit expressions for the classical 
energies of all the plateaux are given in the appendix 
VIII. Additionally we display all the relevant classical 
structures in Fig. 4. 

We obtain a sequence of plateaux at densities 1/9, 
2/15, and 1/6 in the range 0.5 < J'/J < 0.68 using 
Eq. (1) in full agreement with the results given in Ref. 18 
for J'/J — 0.5. The details of the calculation will be 
discussed below. It is remarkable that there are only 
two transitions with two multi-intersectional points (see 
Fig. 3). Note that at these two multi-intersectional points 
many other magnetization plateaux are degenerate which 
we do not display in Fig. 3 for clarity reasons. The energy 
of the 1/8-ca structure shown in Fig. 3 just intersects the 
other curves but the 1/8 plateau is not favored. The two 
plateaux 1/8-diamond and 1/8-tilted are almost degen- 
erate and slightly above the 1/8-ca plateau. The exis- 
tence of the 1/8 plateau in the frustrated quantum mag- 
net SrCu 2 (B0 3 )2 can therefore not be explained solely 
by the effective interactions 18 . The origin must be a con- 
sequence of quantum fluctuations which either originate 
from the Shastry-Sutherland model itself or from addi- 
tional magnetic terms like the DM-interaction. Below, we 
will indeed show that the DM-interaction plays a central 
role for the appearance of the 1/8 plateau with a diamond 
unit cell. 

Although the sequence of plateaux is unchanged when 
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J'/J is increased, the structure of the classical 1/6 
plateau becomes different for J'/J > 0.67. This is a 
direct consequence of the fact that the two-body inter- 
action V 3 ' strongly increases for large J'/J compared to 
the competing terms V4, V5, and V7 (see Figs. 1 and 2). 
It therefore becomes attractive to realize structures con- 
taining the latter interactions and to avoid V 3 '. This is 
the reason why the structure of the 1/6 plateau changes 
for J'/J > 0.67 already in the CA as can be seen in the 
difference between Fig. 3(a) and Fig. 3(b). 

Most interestingly, the phenomenological theory of 
the experimental data 24 proposes the 1/6-stripe plateau 
which does not contain any V 3 interaction but is built by 
V4 and V5 interactions similarly to the 1/6-square and 
the 1/6-new structure shown in Fig. 3(b). All the three 
plateaux have exactly the same classical energy. There- 
fore, we have found strong evidences that the ratio J'/ J 
must be rather large being close to the phase transition 
J'/J 0.67. Let us remark that the accuracy of the 
value J' I J — 0.67, where the 1/6 structure does change 
in the classical limit, of course depends on the accuracy 
of the two-body interactions V$ which might lead to a 
little shift of this value. 



III. QUANTUM FLUCTUATIONS 

In the following we discuss the quantum fluctuations 
which are relevant for the magnetization plateaux at low 
densities. We study first the kinetic terms of the pure 
Shastry-Sutherland model. Afterwards we discuss the 
effect of additional magnetic couplings like DM interac- 
tions or Heisenberg exchanges to more distant neighbors. 



A. Kinetic terms in the Shastry-Sutherland model 

It is known that there exists only one very weak stan- 
dard hopping process t 2 in the Shastry-Sutherland model 
due to the strong frustration 1 ' 27 . This kinetic process 
represents a hopping over the diagonal (see Fig. 1) which 
starts only in order six perturbation theory with a very 
small prefactor 

%^ = * 2 £ £St +a S,. (2) 

a— ±rc±y j 

We have calculated this hopping amplitude up to order 
17 in J'/J using Takahashi's degenerate perturbation 
theory 28,29 . Note that one would get exactly the same 
series using pCUTs. For the given problem, Takahashi's 
expansion is more efficient for one-particle properties and 
one is able to determine a higher perturbative order. The 
extrapolated series are shown in Fig. 5. The amplitude 
is very small in a broad range of couplings J'/J < 0.6. 
Nevertheless, in the most important range J'/J > 0.6 
we observe a rapid increase of the hopping amplitude ti- 
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FIG. 3: Lowest classical energies E C \/(JN) per dimer of all 
plateaux in the density regime 1/9 < n < 1/6 as a function 
of fi/J for (a) J'/J = 0.5 and (b) J'/J = 0.68 using Eq. (10) 
setting all kinetic terms to zero. The energies of the different 
2/15 respectively 1/6 structures shown are exactly degener- 
ate. The circle on a dimer denotes the presence of a particle. 



The coupling is found to be -0.015 < t 2 /J < -0.005 for 
J'/J = 0.65 - 0.70. 

Let us remark that the dominant kinetic terms in the 
effective model are actually so-called correlated hopping 
terms where one triplon is able to hop if another particle 
is present but remains static 4,18,27 . Indeed, correlated 
hopping processes arise already in order two perturba- 
tion theory. But such terms are not important in the 
limit of low densities. Here triplons in the Wigner crys- 
tals are rather far apart and the effects of correlated hop- 
ping is suppressed. As a consequence, the importance of 
correlated hopping increases with increasing density. Be- 
low, we show that correlated hopping is unimportant for 
n < 1/8 but it might be of relevance for n = 1/6. Clearly, 
at larger densities like n = 1/4 correlated hopping is ex- 
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FIG. 5: One-triplon hopping amplitude t2/J as a function of 
J' I J. Solid line represents the bare series of order 17 while the 
dashed curves correspond to various dlogPade extrapolations. 



FIG. 4: Classical plateaux which are relevant at low densities. 
The unit cells of the different structures are shown in dark 
gray (blue). Note that not all plotted structures are realized 
within the CA. This includes all plateaux at density n = 1/8. 
Additionally, for the structures 2/15-b2 (2/15-big) is built 
alternatingly by three 1/8-diamond (three 1/8-tilted) and one 
1/6. Each 1/8-stripe is shaded in light gray. 
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pected to be essential 30 . 



B. Additional magnetic couplings 

The pure Shastry-Sutherland model is a good but 
not perfect microscopic model for the description of 
SrCu2 (B03)2. It is known that a finite DM interaction 
is present which is expected to be a few percent of the 
dominant nearest- neighbor Heisenberg exchange J 9 > 31 ^ 34 . 
Additionally, also Heisenberg exchange interactions be- 
tween next-nearest neighbor dimers might be of a similar 
order 33 . 

It is important to realize that all of these interactions 
are typically unfrustrated and, as a consequence, will 
pump kinetic energy into the system. Here we aim at 
treating the first-order effects on the effective low-energy 
model. This is expected to be a reasonable approxima- 
tion since all the additional couplings are small pertur- 
bations of the order J/100. Let us stress that these per- 
turbations can nevertheless be important at low densities 
because perturbatively the two-body interactions stabi- 



FIG. 6: Illustration of the additional magnetic couplings be- 
yond the ones included in the Shastry-Sutherland model: (a) 
inter-dimer DM interaction D z in z-direction. (b-c) Heisen- 
berg exchanges to next-nearest neighbor dimers are denoted 
by J c (chain-like coupling) and J\ (ladder-like coupling). 



lizing the low-density plateaux appear in a rather high 
order in J'/J as discussed in detail above. 

We start by discussing the effects of DM interactions 
D ■ Si x Sj which turn out to be the most relevant correc- 
tion for SrCu2(BC>3)2. One has to distinguish between 
intra- and inter-dimer DM interactions. 

Taking the DM interaction as a perturbation of the 
bare Shastry-Sutherland model, it can be easily checked 
that the intra-dimer DM interactions leads to a mixing 
of singlet and triplet states on the single dimers. The 
effective low-energy model is therefore only affected in 
second order perturbation theory for this perturbation. 

This is different for the inter-dimer interaction between 
two dimers. We assume that the interaction is largest 
in z-direction and that it has a different sign for the 
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two interactions between nearest-neighbor dimers as il- 
lustrated in Fig. 6 31 ~ 33,35 . In the effective model the 
inter-dimer DM interaction results in a nearest-neighbor 
hopping t\ — D z /2 which reads 



H, 



cff,t, 



E 

y=±y; 



3 



(-1) 



(3) 



where (— 1) J positive (negative) corresponds to effective 
sites j representing a vertical (horizontal) dimer as shown 
in Fig. 1. Additionally, we note that the inter-dimer DM 
interactions D x and D y do not give any boson-conserving 
terms in the effective model in leading order. These cou- 
plings result in processes of the form (b^b^b +h.c). Since 
the -D x .y components are expected to be smaller than D z 
and since they change the terms in the effective model 
only in order D x y ( J' / J) n with n > 1, effects of the trans- 
verse components of the inter-dimer DM interactions are 
expected to be small. 

Another aspect for a microscopic description of the ma- 
terial SrCu2(BC>3)2 are additional Heisenberg exchanges 
to more distant dimers. The couplings to next nearest- 
neighbor dimers are denoted by J c and by J;. They 
are shown in Fig. 6. One can easily see that both cou- 
plings contribute in first order to the already discussed 
diagonal hopping element ti present in the pure Shastry- 
Sutherland model 



t 



total 



J 



Jc 

4 



(4) 



Recently, ab-initio calculations estimated J c /J ps 0.023 
and Ji '/J ps 0.008 for J'/J ps 0.56 33 . The latter ratio 
for J'/J is likely too small. We nevertheless trust the 
order of magnitude for the additional Heisenberg cou- 
plings. Then the total diagonal hopping t^f tal is expected 
to be similar to the amplitude of the pure Shastry- 
Sutherland model. 

Altogether, the most important kinetic terms at low 
densities correspond to a nearest-neighbor hopping t\ 
originating from the inter-dimer DM interaction and to 
a next-nearest neigbor hopping ti already present in the 
pure Shastry-Sutherland model. Both couplings are ex- 
pected to be of the order J/100 for a realistic coupling 
ratio J'/J ps 0.65. 

Let us finally stress that the effective Hamiltonian to 
first order in the additional magnetic interactions has 
still the same symmetries as the pure Shastry-Sutherland 
model in a field. This is especially true for the U(l) sym- 
metry which is only broken once order 2 processes in 
the DM interactions are taken into account. As a conse- 
quence, in our effective low-energy description it is still 
possible that the model displays true superfluid or super- 
solid phases where the U(l) symmetry is spontaneously 
broken at zero temperature. 

In the CA, the energy of the magnetization plateaux 
is independent of the kinetic terms. Thus, the quan- 
tum fluctuations induced by the just discussed kinetic 
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FIG. 7: Illustrations of the regions lZ r having the dimension 
dim(7?. r ) = 25. Left: Figure displays two regions 72. i and 72.2 
of two triplets (black dimers) having an overlap O12 (shaded 
area). Right: Illustration of the regions 1Z T for the classi- 
cal plateau at density n = 1/9. The regions with an odd 
respectively even number are periodically equivalent. Note 
that regions lZ r with r odd and r even display exactly the 
same behaviour due to rotational symmetry of the lattice. 



processes are not at all captured by the CA for the Mott 
insulating phases. In the following we describe our mean- 
field approach aiming at treating the quantum fluctua- 
tions on the different magnetization plateaux. 



IV. APPROACH 

We are interested in studying the effects of the kinetic 
terms t\ and t% on the classical magnetization plateaux. 
Let us stress that we do not want to describe the melting 
of the Wigner crystals clearly present for large kinetic 
terms. We focus on the competition between different 
low-density plateaux because experimentally a sequence 
of plateaux is found at densities 1/8, 2/15, and 1/6. 

The basic idea is that all symmetrically equivalent 
triplons of a Wigner crystal benefit from the kinetic pro- 
cesses in an identical fashion. Additionally, the kinetic 
processes ti and £2 are small compared to the largest two- 
body interactions V\ and V3 and therefore quantum fluc- 
tuations of triplons inside a Wigner crystal are expected 
to be rather local (but not necessarily small) about their 
classical positions. Consequently, we use a mean-field 
approach. 

We split the Hilbert space into periodically equiva- 
lent finite regions 7Z r (see Fig. 7) and we investigate 
the dynamics (quantum fluctuations) of a single triplon 
inside these regions assuming that all other triplons of 
the Wigner crystal remain static. Physically, this corre- 
sponds to an effective one-body problem in a static exter- 
nal potential which is given by the two-body interactions 
V$ of the surrounding triplons. The kinetic hopping am- 
plitudes ti and t% give rise to hoppings of triplons around 
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their classical position leading to an overall energy re- 
duction of the Wigner crystal. This procedure is iterated 
selfconsistently until convergence is reached. 
To be specific, the full effective Hamiltonian 



S j a=±x; j 

a=±y; 

+ t ? £ E^+Ji-(M-Mo)E^ (5) 

a=±x±y j j 

= £ (if W + + ^ - (/, - Mo) E % 

r j 

is split into three different parts. We stress that each 
lattice site is at least part of one region lZ r , i.e. different 
regions are allowed to overlap. The overlap of one region 
lZ r with all other regions lZ r ' is denoted by O r y . For- 
mally, if a dimer j is contained in the overlap j £ O r y , 
it follows j €E lZ r and j € !Z r ' with r' =/= r. 

The explicit expressions of the three terms in the 
Hamiltonian are given by 

+«i E E (-^Wi 

a=±y; j,j+a£O r r > 

+*2 E E Wi (e) 

a=±x±y jeTZr 



E E (-tfsJ+cA 

a=±y\ j+aeO r y 

+^ E E 

a=±i;jgO rr , J 
a=±y; 

+*2 E E S 1+<A 

a=±x±y jg1Z r i 

]+aeO Tr , 

E E .v 'v..V ( 7 ) 



^ = «i E E i-VVS+Jj 

a=±x- j<£R r , 
a=±y; j+a £K r 

+ t * E E (8) 

a=±x±y j^TZ r , 



The first term H ( r ) of the effective Hamiltonian given 
by Eq. (6) contains all hopping processes where the initial 
and the final dimer is only part of region lZ r . Addition- 
ally, the sum over j of the interaction part is restricted 
to the dimers of region lZ r which do not belong to any 
other region. The second part fl" ( e W'5 includes all kinetic 
processes taking place in the overlap O r y of regions lZ r 
and lZ r i . Similarly, this part contains the interactions 
Vsiijnj+s where the index j is a dimer belonging to the 
overlap of regions lZ r and 7Z r > . Let us remark that Nj 
denotes the number of regions lZ r containing dimer j 
which ensures no double counting is done. Finally, the 

~ (r r') 

third part of the Hamiltonian i^in takes into account 
all kinetic hopping terms from region lZ r to region TZ r i 
such that initial and final dimers are not part of the over- 
lap O r . r ' of different regions. 

In the following we neglect the term ^ , because 

we want to study an effective one-body problem, i.e. we 
study the one-triplon quantum fluctuations of one region 
assuming that all other regions remain frozen. Physi- 
cally, this is reasonable because the triplons are expected 
to only fluctuate strongly in the close vicinity of the clas- 
sical positions of the Wigner crystals due to the presence 
of strong repulsive density-density interactions in the ef- 
fective Hamiltonian. 

The next step is to decouple the two-body interaction 
in the Hartree channel 

hjUk ~ ^fij{h k ) - (hj)(h k ) . (9) 

It is important to study how the hardcore constraint is 
violated in the effective mean-field description and how 
one possibly can correct such violations. Clearly, there 
is no conflict in applying Eq. (9) to the first term 
because no overlap is involved. This is different for the 
second term H^° r . r '\ Here it is in principle possible to 
violate the hardcore constraint by placing the particle 
of the region under consideration to a dimer where the 
probability distribution of another particle from a differ- 
ent region is close to one. As a consequence, two particles 
would be on top of each other which is energetically fa- 
vored because the potential term Vs vanishes for 6 = 0. 
This is clearly unphysical. Furthermore, it is also prob- 
lematic that the probability of such a hopping process 
would be equal to one (see Figs. 8(a-b)). 

We are aiming at partially repairing this violation of 
the hardcore constraint via two renormalization factors, 
one for the kinetic part and one for the interaction part. 
The hopping of a particle should depend on whether 
there already is a finite probability for the presence of 
another particle on the involved dimers. Consequently, 
we put the renormalization factor (1 — (n,-))(l — (fij +a ) 

for the hopping term bj +a bj. The second renormal- 
ization factor 1/(1 — (nj)) corrects the potential terms 
Vshj(nj+$) where the dimer j belongs to the overlap 
O r , r ' ■ There are three reasons to add a factor 1/(1— (hj)) 
for the interaction part: i) If two particles are placed 
on the same dimer, the potential should be infinite due 
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FIG. 8: Qualitative illustration of the renormalization fac- 
tors due to the overlap of different regions, (a): We consider 
the initial situation that a particle on dimer hops to the 
dimer 1 which is part of the overlap with a second region. 
All dimers contained in this second region are shaded in dark 
grey (blue) and the number attached on top of these dimers 
correspond to the probablity (density) that the particle in 
this second region is on a specific dimer. (b): If all renor- 
malization factors are absent, the particle which has moved 
to dimer 1 is unaffected by the 10% probability that another 
particle of the second region is already present on dimer 1 . As 
a consequence, the particle does not pay the potential energy 
resulting from this 10% of the density distribution on dimer 1. 
Because the hopping amplitude is iti it would be attractive 
to hop on dimers with a large density of another particle in 
order to avoid high potential barriers. This is clearly unphys- 
ical. (c)-(d): These two figures illustrate the same situtation 
but with the two renormalization factors for the kinetic part 
and for the interaction part of the mean-field Hamiltonian. 
The hopping amplitude it\ for the hopping from dimer to 
dimer 1 is renormalized to (1 — (fi\))it\. At the same time, 
this leads to a change of the density distribution of the par- 
ticle inside the second region. The density expectation value 
(fik) of the particle in the second region is effectively increased 
by <n*)/(l - (m)). 



(a) Density 



(b) Magnetization 



FIG. 9: Schematic illustration of the effective observables 
S^ S a' z and <§2 ■ Left figure represents a local density h from 
a (mean-field) state of the low-energy description. The right 
figure corresponds to the distribution of the local magnetiza- 
tion in the effective model originating from this finite density. 
The radius of the circles is proportional to the square root of 
the density (left) / local magnetization (right) on the dimer 
(left) / spin sites (right). In the right figure, filled (empty) 
circles denote a local magnetization pointing outside (inside) 
the plane. 



As mentionned above, we study the quantum dynam- 
ics of one particle in a region TZ ro assuming that all other 
particles remain frozen, i.e. these particles are not allowed 
to hop and we replace the density operators n by expec- 
tation values (n) for these particles. Consequently, the 
many-body problem is replaced by an effective one-body 
problem in the finite Hilbert space lZ ro . We therefore 
obtain the following mean-field Hamiltonian 



(mf) 



off 



JN r 



n mf tn mf 



(m-mo)E^ ( 10 ) 



where N r corresponds to the total number of regions and 
the sum in the last term is over all dimers j contained 
in region TZ rQ - The explicit expressions for the first two 
terms are given by 



to the hardcore constraint. This is exactly the case 



for 



= 1. ii) If two particles do not overlap, then 



the renormalization factor should be 1 which is true 
for (hj) — 0. iii) If the particle under consideration 
in region lZ r is placed on a dimer j being part of the 
overlap O r y and if, simultaneously, the density of 
the particle in region lZ r r is finite on this dimer, the 
effective potential in the mean-field description should 
increase on this dimer. This is again guaranteed by the 
renormalization factor (see Fig. 8(c)-8(d)). In the follow- 



ing we use Cj = 1 



in order to lighten the equations. 



s jen ro 
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(11) 
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and 



E ^ E 7^ ( ~ 9 ("i>("i+*> + % ) 

+ «1 E E (-^^Ci+aftJ+Jj 
a=±x; j+a£lZ ro 
a=±y; j£O rQ y 

+t2 e E r / r /- ■ (12) 

a=±x±y j+a£lZ ro 

One is therefore left with the following self-consistency 
equations for the mean- field parameters {hj) with 

(^|StS.|^ m )) = (n.)M (13) 

^cff f) (<%> M ) l^ (m) > = e( m )|0( m+1 )) , (14) 

which are solved iteratively. The index m refers to the 
m-th iteration step. Let us stress that the one-body prob- 
lem for the particle in a region is effectively taken place 
in the finite Hilbert space lZ ro because all other degrees 
of freedom are considered to be frozen. Consequently, 
the diagonal elements (j\H^^ \j) with j 6 lZ ra of the 
resulting matrix only depend on the two-body interac- 
tions. Here \j) refers to the one-particle state where 
the particle in region lZ ro is located on dimer j with 
(hj) = 1 while all other particles remain in their (de- 
localized) state of the previous iteration. For the first 
iteration (m = 0), we choose a purely classical state for 
the particles inside the different regions, i.e. the density 
distribution of a particle in one region is {hj) = 1 for one 
specific dimer j. The classical energy E c \ corresponds 
to the energetically lowest configuration of the potential 
terms Vs. Then the triplon sits on its classical position 
jd with (hj ol ) — 1. Thus one obtains the classical energy 

per dimer E cl /N = {N r /N){j c i\H^ f) \j cl ) where N r /N 
just represents the density of the system. 

The non- vanishing off-diagonal matrix elements are ei- 
ther oc ±iti or oc i 2 depending on the dimers which 
are involved. Since all matrices are small, they can be 
solved readily by exact diagonalization. The resulting 
lowest eigenvalue £K rof ) /(JN) corresponds to the mean- 
field energy per dimer of the whole system. The associ- 
ated eigenvector contains the information how all peri- 
odically equivalent triplons delocalize inside the Wigner 
crystal. 

Up to now, we have considered an effective one-body 
problem where every periodically equivalent region dis- 
plays the identical behaviour. We therefore need to cal- 
culate the mean-field parameters respectively the expec- 
tation values of the density operator only once per itera- 
tion. This is surely right for the 1 / 8-diamond plateau (see 
Fig. 4) . But if we consider other plateau structures hav- 
ing more than one triplon per unit cell, we have to take 



care that not every particle is related by symmetry (e.g. 
the structure 2/15-rhomb shown in Fig. 4). We therefore 
have to iterate additionally over all different regions us- 
ing the results of the last iteration for the particles of the 
other regions. 

We therefore obtain the plateau energies and the wave 
function out of the mean-field procedure. We calculate 
the energies and wave functions with ten different start- 
ing configurations. Eight of them are randomly initial- 
ized, one uses the classical solution as input, and one is 
initialized by setting particles in a diagonal fashion into 
the unit cell. The latter is motivated by the observa- 
tion that diagonal configurations often converge rather 
rapidly to the global energy minimum of the mean-field 
equations. The resulting wave function can be used 
to calculate the derealization of the triplons inside the 
Wigner crystals by taking the density expectation value 
of the mean-field ground-state wave function. Let us 
stress that this derealization is a consequence of the 
quantum fluctuations induced by the kinetic processes 
of the effective Hamiltonian. Additionally, one has the 
quantum fluctuations arising from the fact that the effec- 
tive hardcore boson model (except the contribution from 
the DM interaction) has been derived as a high-order se- 
ries expansion from the Shastry-Sutherland model, i.e. 
observables like the density or the local magnetization 
has to be evaluated in the same basis 36 ' 37 as already 
done for the purely classical solution of the effective 
model 18 . Let us mention that the spin-density profile 
has been also calculate for various structures using exact 
diagonalizations 1 9 . 

Here we are especially interested in the local magne- 
tization on each spin S^' z and S^ z of the Shastry- 
Sutherland lattice, because this is the experimentally rel- 
evant quantity for NMR measurements. We therefore 
have calculated the following effective observables 

@i j Z ~ E ( 15 ) 

a 

z - E > ( 16 ) 

a 

as it has already been done in Ref. 18 for the pure 
Shastry-Sutherland model. The coefficients A a and B a 
are determined up to order 10 in J'/J and we have ex- 
trapolated the series using Pade extrapolation. The local 
magnetization is then obtained on each spin-site as illus- 
trated in Fig. 9. Physically, there are two mechanism 
leading to a derealization of the triplons inside the mag- 
netization plateaux. First, the quantum fluctuations of 
the Shastry-Sutherland model itself already present for 
the purely classical Wigner crystals 18 and, second, the 
quantum fluctuations induced by the kinetic processes of 
triplons about their classical positions inside the magne- 
tization plateaux. 

The just presented mean-field theory is expected to 
work well for plateau structures which are not too far 
from classical plateaux of frozen single triplons. But 
there are obvious limitations of our approach which are of 
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(b) Mean-field energies E m f/(JN) per dimer 



FIG. 10: (a): The energy gain per dimer E &a x n /{JN) is 
plotted as a function of ti/J and t%/ J for density 1/8 at 
J'/J = 0.65. The three considered structures are 1/8- 
diamond (triangle down), 1/8-tilted (circle), and 1/8-ca (tri- 
angle up) . The combination of both kinetics terms favors the 
structure 1/8-diamond with a diamond unit cell. For other 
J'/J values (not shown), the energy gain E ga .i n /(JN) has a 
very similar behaviour, (b): The mean-field energy per dimer 
Emf/(J N) is plotted for 1/8-diamond (triangle down), 1/8- 
tilted (circle), and 1/8-ca (triangle up) as a function of ti/J 
and t2/J. Note that in the shown parameter regime always 
one of these three structures corresponds to the global mini- 
mum of the mean-field energy. 



two different types: a) The kinetic terms are effectively 
too large compared to the involved repulsive interactions. 
Then a superfluid (supersolid) solution is expected and 
any crystalline solution is unphysical. Clearly, this has to 
happen in the dilute low-density limit. In our approach 
targetting Mott insulating phases, the involved repulsive 
interactions typically stabilizing the Wigner crystals are 



negligibly small because the distance between triplons be- 
comes large at low densities. As a consequence, already 
small values of the kinetic terms are expected to melt 
any formation of a Wigner crystal. In the mean-field cal- 
culation, one then observes a negative mean-field energy 
Emf < at fi = fiQ which we consider as an indicator for 
the presence of superfluid phases. The mean-field solu- 
tions having E m f < are unphysical and we will denote 
such solutions in the following to be of a-type. f3) Our 
mean-field treatment is expected to breakdown when the 
reduction to an effective one-body problem fails. Later, 
we observe in certain situations that the mean-field so- 
lution corresponds to a Wigner crystal where triplons 
inside their region completely delocalize among various 
dimers. Such a behaviour is likely an artefact of our 
mean-field solution and it is therefore unphysical. In any 
case, one cannot expect that such strongly delocalized 
plateaux are well described on the mean-field level. We 
will denote such problematic solutions in the following 
to be of /3-type. To be specific, we classify a mean-field 
solution as /3 if the second-largest density of a triplon on 
a dimer inside a region exceeds the value 0.25. 



V. RESULTS 

In this section we present the results obtained by the 
mean-field calculation. The numerical setup we use is to 
iterate over all unit cells up to 100 dimers getting 257 
different plateau densities for n < 1/6. We focus on 
the range 0.6 < J'/J < 0.7 with < t\jJ < 0.02 and 
< —til J < 0.02, which is expected to be the relevant 
regime for SrCu2(B03)2- 

In the following, we study first the favored plateau 
structures for a fixed density. This is done for n = 1/8, 
n = 1/6, n = 2/15, and n = 1/9. Afterwards we dis- 
cuss the full low-density phase diagram of the extended 
Shastry-Sutherland model in Eq. 5. For all cases we give 
a detailed comparison to the physics of the frustrated 
quantum magnet SrCu2(B03)2- 



A. 1/8 plateau 

A triplon can well delocalize to neighboring dimers if 
the potential difference is small. The largest repulsive 
two-body interactions are V\ and V3. All dimer sites 
which involve these two interactions represent therefore 
a sizable energy barrier and fluctuations to these dimers 
are well suppressed. To understand which kind of 1/8 
plateau is favored, let us start by investigating the differ- 
ent structures at fixed density n — 1/8 in our numerical 
setup. To be concrete, we compare the structures 1/8- 
diamond, 1/8-tilted, and 1/8-ca (see Fig. 4) as well as 
possibly other structures in case they have the lowest 
energy. 

We consider first the energy gain .Egain for these three 
structures defined by the difference between the classi- 
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(c) 1/8-ca 



FIG. 11: Density distribution of the mean-field solution for 
density 1/8 at J'/J = 0.65, h/J = 0.01, and t 2 /J = -0.005. 
Note that the densities of all dimers are rounded to the third 
decimal digit. The suppression of certain fluctuation channels 
are clearly visible. The radius of the plotted circles is propor- 
tional to the square root of the density on each dimer. The 
red circles denote densities which are larger then 0.5. 



cal energy E c \ and the converged mean-field energy E m { 
including quantum fluctuations 



E ■ 1 



(17) 



It is remarkable that the 1/8-diamond structure has in 
most cases the highest energy gain (see Fig. 10(a)). Qual- 
itatively, we find that increasing t\/J gives basically the 
same high energy gain for the 1/8-diamond and for the 
1/8-tilted structure. In contrast, the energy gain for the 
classical plateau structure 1/8-ca is rather small. Fur- 
thermore, increasing i 2 /J results in an energy splitting 
between the preferred plateau with a diamond unit cell 
1/8-diamond and the other two structures at density 
n = 1/8. 

These findings are well understood by looking at the 
potential differences between the classical position of the 
triplons inside the Wigner crystals and the dimers which 
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FIG. 12: Minima of the mean-field energy are plotted as a 
function of t\/ J and of t 2 / J. The dashed blue line separates 
the parameter regime where the classical structure 1/8-ca is 
favored from the regime where the 1/8-diamond structure is 
realized. This line is obtained by fitting the lowest energies of 
1/8-diamond and 1/8-ca along the displayed grid. Symbols a 
(/3) represent mean-field solutions of a- type (,3-type). Mean- 
field solutions which differ by less than 10 _4 J are considered 
to be degenerate. In this case all such structures are displayed. 



can be reached via one hopping t\ or one hopping £ 2 (see 
Fig. 11). For a finite t\ and = 0, only the two plateaux 
1/8-diamond and 1/8-tilted are such that all four neigh- 
boring dimers do not involve the large repulsive interac- 
tions Vi or V3 (compare Figs, ll(a-b)). In contrast, the 
classical structure 1/8-ca is stabilized by V 3 ' interactions, 
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and therefore two out of the four nearest neighbors in- 
volve the largest repulsive interaction V\ and quantum 
fluctuations to these dimers are suppressed. One then 
expects that a finite t\ favors the 1/8 structures with 
diamond and tilted unit cells in a similar fashion. This 
is different for the other limiting situation, i.e. a finite 
ti and a vanishing t\. Here it is only the 1/8-diamond 
structure where all four (equal) fluctuation channels do 
not involve large interactions terms. In contrast, for the 
other two structures minimally one out of four fluctua- 
tion channels are suppressed due to the interaction V 3 . 
This explains the energy splitting between the tilted and 
the diamond structure triggered by i 2 as discussed above 
(see Fig. 10(a)). 

Next we compare the ground-state energies E m { of 
the different plateau structures at density 1/8 which are 
shown for J'/J = 0.65 in Fig. 10(b). Interestingly, it is in 
principle only the 1/8 structure with diamond unit cell 
which is stabilized besides the already expected classi- 
cally realized 1/8-ca structure. The unphysical structure 
denoted by f3 for J' / J — 0.68 is of /?-type, i.e. triplons are 
delocalized completely on two different dimers. Further- 
more, it is the hopping to nearest-neighbors t\ originating 
from the DM interaction which is the driving force for the 
1/8-diamond structure found within the phenomcnolog- 
ical theory for SrCu 2 (B0 3 )2 24 . Indeed, we checked ex- 
plicitly that the next-nearest neighbor hopping i 2 alone 
is not sufficient to stabilize the structure with a diamond 
unit cell. These trends are present for all ratios of J' /J 
as long as J'/J is large enough (see Fig. 12). Indeed, 
even for a sizable J'/J = 0.6 only a tiny region of 1/8- 
diamond is found. Instead, we observe that in most cases 
mean-field solutions of a-type have the lowest mean-field 
energy. 

As stated above, solutions of a-type indicate a melting 
of the Wigner crystals which is expected once the kinetic 
terms are large enough to win over the potential terms. In 
our mean-field treatment focusing on crystalline phases, 
an indication for such a melting transition to superfluid 
phases is a negative mean-field energy E m { at /i = fio as 
can be seen for example in Fig. 10(b) for ti/J = 0.02 
and tij J = —0.02. The melting transitions of the dif- 
ferent Wigner crystals have to occur for smaller values 
of ti and ti when the ratio of J'/J is reduced, because 
then also the repulsive interactions stabilizing the mag- 
netization plateaux are smaller. This is exactly what we 
observe in the behaviour of the mean-field energy which 
becomes negative at /i = for smaller values of the ki- 
netic terms if J'/J is reduced (see Fig. 12). The presence 
of a 1/8 plateaux with diamond unit cell in the frustrated 
quantum magnet SrCu 2 (B03) 2 points therefore to a ra- 
tio J'/J rj 0.65 very close to the phase transition point 
in the Shastry-Sutherland model. The latter finding is 
further confirmed in the next subsection when studying 
the magnetization plateau at fixed density 1/6. 

Finally, we want to see whether the above discussed 
1/8 structures (especially the one with a diamond unit 
cell) are indeed realized in the phase diagram. To do so 
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FIG. 13: The figure displays for J'/J = 0.68 the range in 
the chemical potential A/x/J where a 1/8 plateau is realized 
in the full low-density phase diagram compared to all other 
plateau structures at different densities. The stabilized 1/8- 
plateau changes its structure from the classical 1/8-ca struc- 
ture (dashed lines, present for small values of ti / J and — t-z / J) 
to the l/8diamond structure (shown as solid black line) when 
both hoppings t\/J and —tij J are of the order 0.01. 



we compare all different plateaux at different densities as 
a function of the chemical potential fi/J and we plot in 
Fig. 13 the ^-range 

A/i Hcnd — Mstart „n 

T = 1 (18) 

where a 1/8-plateau is realized as a function of t\j J and 
t 2 /J for J'/J = 0.68. Here /i on d (/^start) denotes the 
end (beginning) of the plateau at density 1/8 in the low- 
density phase diagram. The quantity Afi/J corresponds 
therefore to the width of the plateau. Let us mention 
again that the 1/8 plateau is not favored in the classi- 
cal limit, i.e. for t\ = t% = 0. It is therefore remark- 
able that already for small kinetic hopping terms a 1/8 
plateau with the 1 /8-ca structure is present in the mean- 
field phase diagram. But one should keep in mind that 
the transition point between the plateau at 1/9 and at 
2/15 in the classical limit is highly degenerate, i.e. many 
structures including the 1/8-ca have the same classical 
enegy for this specific chemical potential. Thus small ki- 
netic terms are able to lift this degeneracy and to favor 
the 1/8-ca plateau in a tiny /i- range depending on the 
size of the kinetic terms. 

More importantly, we find that the 1/8 plateau with di- 
amond unit cell is present in the phase diagram for a wide 
range of hopping parameters including the expected val- 
ues for SrCu 2 (B0 3 ) 2 : h/J « 0.015* and t 2 /J « -0.01. 
It is indeed the combination of both kinetic terms which 
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is responsible for the stabilization of the 1/8-diamond 
plateau. 

Altogether, our mean-held results for density 1/8 are in 
very good agreement with the experimental findings and 
the phenomenological theory for SrCu2(B03)2 24 . Quan- 
tum fluctuations induced by the kinetic hopping terms t\ 
and ti are essential to obtain a sizable 1 /8 plateau with a 
diamond unit cell. Furthermore, the coupling ratio J'/J 
must be rather large J'/J ~ 0.65 in order to prevent 
a melting of the Wigner crystal due to the kinetic pro- 
cesses. 



B. 1/6 plateau 

The other prominent low-density plateau observed ex- 
perimentally in the compound SrCu2 (1303)2 is at density 
1/6 which we focus on in this subsection. Interestingly, 
the phenomenological interpretation of the NMR data 
on SrCu2 (1303)2 24 yields a structure for the 1/6 plateau 
which is in disagreement with the plateau 1 /6-ca found in 
the classical limit for J'/J — 0.5 18 . The latter structure 
is stabilized by V3 interactions. This has to be contrasted 
with the phenomenologically deduced 1/6 plateau which 
is built out of V4 and V5. As stated above, we find ex- 
actly such structures (1/6-stripe, 1/6-square, and 1/6- 
new) naturally already in the classical limit for ratios 
J'/J > 0.67. So again, our results indicate a rather large 
ratio J'/J in SrCu 2 (B0 3 ) 2 . 

The three plateau structures 1/6-stripe, 1/6-square, 
and 1/6-new have exactly the same classical energy. In 
the following we want to study how this degeneracy is 
lifted in our mean-field theory if the kinetic processes t\ 
and i 2 are turned on. To this end we discuss first the 
dominant fluctuation channels for these plateaux (plus 
the structure 1/6-ca) in an analogue fashion as we did 
for the plateaux at density 1/8. One therefore has to 
check for each particle on each structure whether the po- 
tential barrier on the eight dimers which can be reached 
from the classical position by one t\- or one t 2 -hopping 
is large or not. 

In the 1 / 6-ca structure all particles have the same en- 
vironment. For each particle, there are only two dimers 
which can be reached via one ii-hopping from the classi- 
cal position that do not involve the large V\ or V 3 interac- 
tions. All other six fluctuation channels are suppressed. 
This is different for the other structures. The two struc- 
tures 1 / 6-stripe and 1 /6-square stay degenerate when the 
quantum fluctuations of the particles are restricted to 
the eight dimers around the classical position. For both 
structures one finds that four out of the eight fluctuation 
channels do not involve Vi or V3. To be concrete, for 
each particle inside the Wigner crystal there are two out 
of four dimers (for t\ and i 2 ) where the largest poten- 
tial term to be paid is V<z- Finally, the 1/6-new plateau 
has very similar quantum fluctuations as the latter two 
structures. The key difference is the existence of one ad- 
ditional fluctuation channel for the triplons having their 
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FIG. 14: The energy gain per dimer E S3b i n /(JN) is plotted 
as a function of t\/J and t^/J for fixed density n = 1/6 at 
J'/J = 0.65. In contrast to the structure 1/6-new (squares), 
the classical structure 1/6-ca (triangles down) has a lower 
energy gain for all considered parameter values. The structure 
1/6-square (1/6-stripe) is denoted by triangles up (circles). If 
the kinetic terms become large ti/J ~ 0.02 and til J ~ —0.02, 
an unphysical plateau of /3-type has the lowest energy. 



classical location on the vertical dimers inside the Wigner 
crystal (see Fig. 4). This additional channel is reached 
via one t\ hopping to a dimer where one has to pay only 
V4 interactions. The structure 1/6-new is therefore ex- 
pected to gain most from quantum fluctuations. This is 
exactly what we get (see Fig. 14). We therefore find that 
the 1/6-new plateau is favored for most combinations of 
kinetic hopping terms t\ and t2 for J'/J = 0.68. This is 
only different for small values of ti where the structures 
1/6-square and 1/6-stripe have a slightly smaller energy 
(see Fig. 14). 

Let us stress that for too large values of the kinetic 
hopping terms a different Wigner crystal has the lowest 
mean-field energy. As already mentionned above, this 
novel structure is likely an unphysical artefact of our 
mean-field treatment. To be specific, we find a mean-field 
solution where the triplons are not anymore well localized 
at their classical position. In contrast, the particles get 
completely dclocalized among different dimers in their 
regions which is not the behaviour where our mean-field 
treatment is expected to work well. Such mean-field so- 
lutions are of /3-type as introduced in Sect. IV. 

Let us compare our mean-held results with the findings 
deduced from the NMR data of SrCu 2 (B0 3 )2 24 . The lat- 
ter militates in favor of the 1/6-stripe structure. Inter- 
estingly, we find that this plateau is indeed among the 
ones with the lowest energy as long as the ratio J'/J is 
sufficiently large. In this respect our results are in bet- 
ter agreement with the experimental data as the purely 
classical results for J'/J — 0.5 18 . But the structure 1/6- 
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FIG. 15: The structures at fixed density n — 1/6 having 
the lowest mean-field energy E m t are plotted as a function of 
ti/J and t 2 /J for (a) J'/J= 0.65 and for (b) J'/J = 0.68. 
The most important difference between J'/J = 0.65 and J'/J 
= 0.68 is the change from the 1/6-ca structure to the 1/6- 
new structure. Symbols j3 represent mean-field solutions of fi- 
type. Mean-field solutions which differ by less than 10 _4 J are 
considered to be degenerate. In this case all such structures 
are displayed. 



new has a slightly lower energy which one can understand 
due to the larger number of fluctuation channels as dis- 
cussed above. This mismatch between microscopic and 
phenomenological theory originates most likely from sub- 
leading effects not contained in our calculation. Indeed, 
the experimental data do not find evidences for rather 
larger unit cells like for example the structure 1/6-new. 
One possible explanation of this discrepancy is the rele- 
vance of three-dimensional couplings or the renormaliza- 
tion of the two-body interactions V$ due to the additional 
couplings. 

Nevertheless, our mean-field results as discussed so far 
for densities 1/8 and 1/6 are quite promising. The find- 
ings for 1 /8 are in full agreement with experiments and its 
phenomenological interpretation. Additionally, the anal- 
ysis for both densities point to a large ratio J' / J > 0.67 
and one clearly sees the relevance of the quantum fluctu- 
ations we have introduced in our microscopic calculation. 



C. 2/15 plateau 

The NMR data on SrCu2 (1303)2 reveals the existence 
of a third low-density plateau at density 2/15 as predicted 
theoretically for J'/J = 0.5 by the classical solution of 
the effective hardcore boson model 18 . But as for the 
plateau at 1/6, the classical structures 2/15-rhomb and 
2/15-rect are inconsistent with the experimental data. 
Interestingly, the phenomenological theory point to the 
stripe structure 2/15-big 24 , i.e. this Wigner crystal is 
built by a pattern of three stripes of the structure 1 /8- 
tilted and one stripe of the structure 1 / 6-stripe. Let us 
mention that the phenomenological interpretation of the 
NMR data for the 2/15 plateau is very complicated, be- 
cause the experimental signal is very complex due to the 
large unit cell. 

In the following we compare the results of our mi- 
croscopic theory with the above predictions. Most im- 
portantly, our mean-field theory naturally gives a stripe 
structure for the density 2/15. The microscopic reason is 
rather simple. As we have seen in the last sections, the 
1 /8 structure with a diamond unit cell is one of the very 
rare Wigner crystals where all eight fluctuation channels 
to dimers reached by one t\ or one t 2 hopping do not 
involve the large repulsive interactions V\ and V3. As a 
consequence, it seems preferable to construct stripe struc- 
tures above 1/8 which contain the 1/8-diamond Wigner 
crystal as a substructure. One relevant example for such 
a stripe structure with density 2/15 is the structure 2/15- 
b2 illustrated in Fig. 4. One can see that the 2/15-b2 
plateau is formed by a subpart of three 1/8-diamond 
stripes (gray shaded) and one 1/6 subpart which is built 
by V4 and V5 potential terms. It is clear that such stripe 
structures will mainly benefit from the 1/8 component 
since the fluctuation channels for the thin 1/6 stripe are 
mainly suppressed. 

After this discussion, it is obvious that it is less advan- 
tegeous to form stripes by taking 1/8-tilted substructures 
instead of 1/8-diamond for all combinations of J'/J, 
ti/J, and ti/ J . The 2/15 striped crystal formed by 1/8- 
diamond substructures has always a lower mean-field en- 
ergy than the one consisting of tilted substructures. We 
therefore find that the striped 2/15 structure contain- 
ing 1/8-diamond substructures are favored nearly in the 
same region as the 1 /8-diamond crystal itself (see Figs. 12 
and 16). 



D. 



1/9 plateau 



Before we present the resulting phase diagram for real- 
istic values of the kinetic hopping terms, let us first dis- 
cuss our results for fixed density n — 1/9. Neither Ref. 7 
nor Ref. 24 find experimental indications for a plateau 
at density 1/9. In contrast, Sebastian et al have inter- 
preted their high-field torque measurements as evidence 
for a plateau at 1 /9 8 which is clearly in conflict with the 
other measurements. 
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FIG. 16: The structures at fixed density n = 2/15 having the 
lowest mean-field energy E m f are plotted as a function of ti/J 
and t 2 /J for (a) J'/J= 0.65 and for (b) J'/J = 0.68. The 
dashed blue line separates the region favoring the structure 
2/15-rect present for small values of the kinetic terms from 
the regime where the structure 2/15-b2 is realized. This line is 
obtained by fitting the lowest energies of 2/15-rect and 2/15- 
b2 along the displayed grid. Symbols a (/3) represent mean- 
field solutions of a-type (/3-type). Mean-field solutions which 
differ by less than 10 _4 J are considered to be degenerate. In 
this case all such structures are displayed. 



On the theoretical side, the CA 18 finds at J'/J=0.5 
a stable 1/9 plateau and not the one with density 1/8. 
Furthermore, we find a 1/9 plateau with this classical 
structure (see Fig. 4) for a wide range of kinetic couplings 
for all investigated J'/J. One nevertheless expects that 
this 1/9 plateau should melt already for rather small val- 
ues of the induced quantum fluctuations triggered by the 
kinetic hopping processes due to the fact that the 1/9 
plateau is stabilized by the weak repulsive interaction Vq 
appearing only in order 8 perturbation theory. 

An upper bound for this melting transition can be ob- 
tained in our mean-field calculation by determining the 
sign of the mean-field energy at /i = fiQ. This corresponds 
to mean-field solutions of a-type which is shown in Fig. 17 
for J / J = 0.65 and J'/J = 0.68. One clearly sees that 
the 1/9 plateau is about to melt for realistic values of 
the kinetic hopping processes t\ and £2- It is therefore 



FIG. 17: The structures at fixed density n = 1/9 having the 
lowest mean-field energy E m { are plotted as a function of ti/J 
and ti/J for (a) J'/J= 0.65 and for (b) J'/J = 0.68. The 
dashed blue line separates the parameter regime where the 
1/9-plateau has a positive mean-field energy at n — fio from 
the parameter region where the mean-field energy becomes 
negative at /1 = no which is clearly unphysical. This unphys- 
ical behaviour is denoted by a. 

likely that the 1/9 plateau is not a stable Wigner crys- 
tal for realistic parameters because quantum fluctuations 
destroy the Mott insulator. But it is beyond the validity 
of our mean-field approach to answer this question in a 
quantitative fashion. 



E. Phase diagram 

In the last paragraphs we have concentrated on com- 
paring different Wigner crystals at fixed densities 1/8, 
1/6, 2/15, and 1/9 which is motivated by the recent ex- 
perimental NMR data on SrCu2 (1303)2 and their phe- 
nomenological interpretation 24 . Our microscopic mean- 
field theory of the effective hardcore boson model plus ad- 
ditional kinetic terms confirms many aspects of the phe- 
nomenological theory. Altogether, all our results point to 
the following coupling ratios for a microscopic descrip- 
tion of SrCu 2 (B0 3 ) 2 : J'/J > 0.67, h/J w 0.015, and 
til J w —0.01. Interestingly, the magnitude of the ki- 
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netic terms are in very good agreement with estimates 
for the pure Shastry-Sutherland model (see discussion 
on t 2 above) as well as estimates for the DM interaction 
D z /J = 2ti pa 0.02 31 - 35 . The latter is estimated for 
J'/J « 0.635. Obviously, such estimates give the correct 
order of magnitude but do not pinpoint the value quanti- 
tatively. Additionally, our findings of a rather large ratio 
J'/J > 0.67 might result in slightly renormalized fitting 
values. In this subsection we study the full low-density 
phase diagram in this most realistic coupling regime and 
we compare our results to the experimental findings. 

We start by discussing the phase diagram of the pure 
Shastry-Sutherland model which contains intrinsically 
the diagonal hopping term t 2 . It is therefore already 
interesting to compare our mean-field results treating 
quantum fluctuations with the CA used in Ref. 18 giving 
a sequence of magnetization plateaux at densities 1/9, 
2/15, and 1/6 for J'/J = 0.5. 

The density curves for J'/J G {0.6,0.65,0.68} and 
different values for t 2 /J are shown in Fig. 18. One 
clearly sees that the density displays basically the same 
sequence 1/9, 2/15, and 1/6 as the CA 18 . The only ma- 
jor change is for the 1/6 plateau whose structure changes 
for J'/J > 0.67 as discussed above which is consistent 
with the phenomenological theory . Furthermore, it is 
not surprising that already small kinetic terms lift the 
degeneracy of the multi-intersectional points present in 
the classical phase diagram (see Fig. 3). Most impor- 
tantly, the diagonal hopping t 2 alone is not responsible 
for the experimentally stabilized 1/8 plateau. These find- 
ings are well seen in Fig. 18: other plateaux distinct from 
the classical structures at densities 1/9, 2/15, and 1/6 are 
basically not realized. 

Next, we introduce a finite value for the nearest- 
neighbor hopping t\/J — 0.01 and we study again the 
low-density phase diagram shown in Fig. 19. Still, the se- 
quence of plateaux is close to the one already observed in 
the classical limit for J'/J < 0.65. The first real change 
happens for J' / J = 0.68. Here the 1/8 plateau is realized 
with an appropriate width. Nevertheless the structure of 
the 1/8 plateau is still incompatible with the one seen 
in experiments. This suggests that one should consider 
slightly larger values for the nearest-neighbor hopping t\ 
in our mean-field calculation. 

We therefore checked first whether a (realistic) pa- 
rameter range exists for t\ and t 2 where a 1/8 plateau 
with diamond unit cell and a 2/15 plateau with a striped 
structure is realized simultaneously in the phase diagram, 
i.e. where both plateaux are present in the density curve. 
This is shown in Fig. 20. We find that the nearest- 
neighbor hopping ti should be t\/J ~ 0.015. Let us 
stress that it is well possible that such absolute values 
slightly shift due to uncertainties in the extrapolation of 
the potential terms Vs. But it is already encouraging 
that the value for ti deduced from our mean-field theory 
is so close to the one proposed in the literature 9 ' 31-34 . 

Summarizing the above results, we should focus on a 
ratio J'/J > 0.67 and we should consider kinetic cou- 
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FIG. 18: The density n is plotted as a function of /x/J for a 
vanishing nearest-neighbor hopping t\ — for different values 
of J'/J and t 2 /J. The value of the diagonal hopping t 2 is 
consistent with the extrapolated order 17 series of the pure 
Shastry-Sutherland model. Additionally, the most important 
structures are shown. 
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FIG. 19: The density n is plotted as a function of fi/J for 
a nearest-neighbor hopping t\ — 0.01 and different values 
of J' I J and t 2 j 'J. The value of the diagonal hopping t 2 is 
consistent with the extrapolated order 17 series of the pure 
Shastry-Sutherland model for J'/ J > 0.65. Additionally, the 
relevant structures are shown. 
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FIG. 20: Structures with the lowest mean-field energy are 
plotted in comparison to all other densities as a function 
of ti/J and t 2 /J for the density (a) 1/8 and for the density 
(b)2/15. The black symbol x marks a parameter regime where 
both plateaux are realized in the phase diagram. In the rel- 
evant regime ti/J w 0.015 and t 2 /J w -0.01, the 2/15b 2 
plateau is only realized in a small region in contrast to the 
1/8-diamond. Symbols a represent mean- field solutions 
of a-type (/3-type). Mean-field solutions which differ by less 
than 10 _4 J are considered to be degenerate. In this case all 
such structures are displayed. 



plings ti/J 0.015 and t 2 /J~ —0.01. Consequently, 
we plot the density as a function of fi/ J for J'/ J = 0.68, 
h/ J = 0.016 and t 2 /J = -0.01 in Fig. 21(a). Addi- 
tionally, we present our results for the density pattern 
(Fig. 21(a)b-d-f) and for the magnetization (Fig. 21(a)c- 
e-g). The latter is calculated according to Eqs. 15 (see 
also Fig. 9). Most prominently, a 1/8 plateau with dia- 
mond unit cell as well as the 1 / 6- new plateau is realized 
in the phase diagram. The just discussed 2/15 striped 
plateau is also realized but its width is very tiny. 

Furthermore, there are additionally three interesting 
density regimes which one should discuss: i) densities 
below 1/9, ii) densities 1/8 < n < 2/15, and hi) densities 
2/15 < n < 1/6. We start the discussion with the regimes 
i) and hi) which are denoted by a and f3 in Fig. 21(a). 
We believe that our mean-held theory gives unphysical 
results for both density regions, but due to different phys- 
ical reasons. In the low-density region a, we hnd a neg- 
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FIG. 21: (a) Density n is plotted as a function of /i/J for 
J'/J = 0.68, h/J = 0.016, and t 2 /J = -0.01. Plateaux 
shown in red correspond to valid Wigner crystals. Mean-field 
solutions in the regimes a and j3 are unphysical: In a the 
mean-field energies at [i — /10 are negative. In /3 particles 
are completely delocalized among various dimers inside the 
crystalline solution. The structure and the magnetization of 
the plateau at 1/9 are given in (b)-(c). The corresponding 
figure for 1/8 and 1/6 are shown in (d)-(e) and (f)-(g). 



FIG. 22: Figure illustrates the local density of two unphysical 
mean-field states classified as /3 at densities 1/7 (a) and 3/20 
(b) shown in Fig. 21(a). 



ative mean-field energy E m f at /1 = (Xq for all plateau 
structures below 1/9. As explained above, it is very likely 
that these magnetization plateaux melt due to the intro- 
duced quantum fluctuations if one allows also superfluid 
mean-field solutions (or treat the problem beyond mean- 
field theory). In contrast, for the regime j3 between 2/15 
and 1/6, all plateau structures are such that at least one 
particle is completely delocalized between at least two 
different dimers. Two examples of such structures for 
densities 1/7 and 3/20 are displayed in Fig. 22. These 
are therefore clearly structures which are extremely far 
away from perturbed classical plateau structures where 
our mean-field theory is expected to work well. In total, 
we do not trust our mean-field results in the regions a 
and /?. 

The second regime ii) which includes the above dis- 
cussed striped 2/15 structure is different. In fact, the 
same microscopic mechanism giving rise to the 2/15 
plateau, applies to an (in principle) infinite number of 
striped plateaux above 1/8. The specific stripe struc- 
tures together with their magnetizations are shown in 
Fig. 23. Most importantly, our findings agree well with 
the experimental data of SrCu2(BC>3)2 9 ' 24 : One finds 
that the translational symmetry stays broken above the 
1/8 plateau. 



F. Correlated hopping 

Let us finally discuss to what extent so-called corre- 
lated hopping terms can alter the above findings. We 
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have computed all correlated hopping terms starting at 
most in order 6 in J'/ J up to order 14 in J' /J. As 
already mentionned, such couplings allow a hopping of 
a particle when another particle is present and remains 
static. The quantum fluctuations induced by correlated 
hopping depends therefore on the density. The higher 
the density, the more important is the effect of correlated 
hopping for the formation of plateau structures. Let us 
remark that correlated hopping has been also identified 
as a possible driving force for the realization of so-called 
pair superfluids 38 ' 39 . If realized, such phases are expected 
to be present at very low densities in the density range 
which we have classified as a. The description of such 
phases goes well beyond our mean-field considerations. 

Formally, a correlated hopping process is written as 
^ij tMPfrk' i ,e ' a P ar ticle hops from dimer j to dimer i 
with an amplitude t\ ^ k if another particle is present on 
dimer k. The importance of such a correlated hopping 
term for a specific Wigner crystal depends on one hand 
on the strength t\ j k and on the other hand on the po- 
tential barrier on dimer i on which the particle hops. A 
correlated hopping process induces therefore large quan- 
tum fluctuations whenever its amplitude is large and the 
involved repulsive interactions in intermediate states are 
small. 

Let us discuss the most important correlated hopping 
terms for the relevant structures at densities 1/9, 1/8, 
and 1/6 as discussed above. For n = 1/9, there is no 
correlated hopping process starting up to order 6 which 
induce quantum fluctuations on the classical 1/9-plateau. 
One can therefore conclude that correlated hopping has 
no real impact on the mean-field energy of this plateau. 

For the relevant structures at density n — 1/8, there 
exist correlated hopping terms starting in order six 
perturbation theory having an amplitude of the order 
J/1000 in the regime J'/ J « 0.65 (see Fig. 24). In our 
mean- field treatment, these terms renormalize slightly 
the values for t\ and Furthermore, the induced quan- 
tum fluctuations are very similar for the structures 1/8- 
diamond, 1/8-tilted, and 1/8-ca. Altogether, all our find- 
ings for density 1/8 are expected to be unchanged when 
correlated hopping terms are included. 

The situation is different for the structures at den- 
sity n — 1/6. Here correlated hopping processes exist 
which start already in order four perturbation theory. 
The corresponding amplitudes are of the order J/100 for 
J'/ J rj 0.65 (see Fig. 24). The most relevant hopping 
processes are those which do not contain any of the large 
repulsive interactions V\ or V3 in intermediate states be- 
cause then in full analogy to the discussion above fluc- 
tuations from the classical position are not strongly sup- 
pressed. Interestingly, the induced quantum fluctuations 
due to these correlated hopping terms are again almost 
identical for the structures 1/6-stripe, 1/6-square, and 
1 /6-new. We therefore expect that our results for density 
1/6 also hold when correlated hopping terms are taken 
into account. 

The most subtle density regime is the one between 1 /8 
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(b) Magnetization of the stripe structures 



FIG. 23: Figure shows all stripe structures present in the 
low-density phase diagram for J' /J — 0.68, tl/J = 0.016, 
and til J = -0.01 as displayed in Fig. 21(a). All Wigner 
crystals contain stripes build by the 1/8-diamond structure 
(shaded in grey). 



and 1/6 where our mean-field calculation predicts an in- 
finite sequence of stripe structures. Again, there exist 
correlated hopping terms starting in order four perturba- 
tion theory which are of similar magnitude as the kinetic 
processes t\ and t 2 . Consequently, one expects that cor- 
related hopping may have some impact on this hierarchy 
of plateau structures. This might be especially true for 
the density regime just below 1/6 where our mean-field 
calculation stabilizes unphysical states classified as /3 in 
the phase diagram. It would be therefore highly inter- 
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esting to study this part of the phase diagram in more 
detail, because one might expect supersolid phases due 
to the increasingly important correlated hopping terms 5 . 



VI. SUMMARY 

The main motivation of this work is the discrepancy 
between experiment and theory on the understanding 
of the magnetization curve of the frustrated quantum 
magnet SrCu2 (603)2- The most prominent magneti- 
zation plateau at M = 1/8 is currently not under- 
stood microscopically. Instead a sequence of plateaux at 
M = 1/9, M = 2/15, and M — 1/6 is found at low 
magnetization 18 . Interestingly, additional low-density 
plateaux at M = 2/15 and M — 1/6 have been recently 
observed experimentally 24 but their structures disagree 
with theory. 

Microscopically, it is believed that the material 
SrCu2 (B03)2 is well described by the Shastry-Sutherland 
model. In Rcf. 18, an appropriate effective low-energy 
theory for the Shastry-Sutherland model in an external 
magnetic field is derived which is then solved in the clas- 
sical limit. The classical approach is believed to work well 
because quantum fluctuations are strongly suppressed 
due to the strong frustration in the Shastry-Sutherland 
model. The latter is certainly true as long as one is not 
too close to the phase transition point J' / J « 0.7 or 
as long as additional magnetic couplings are negligible. 
Additionally, one assumes that bound states of triplons 
are not essential for the formation of plateaux at low 
densities 25 . 

Here we have studied the effect of quantum fluctu- 
ations on the sequence of magnetization plateaux at 
low densities n < 1/6 which originate either from ad- 
ditional magnetic couplings like the DM-interaction or 
from the one-particle hopping already contained in the 
Shastry-Sutherland model. Both couplings are expected 
to be of the order J/100 at strong coupling J'/ J ps 0.65 
which is the relevant parameter regime for SrCu2(B03)2- 
Physically, the plateaux are formed by freezing hardcore 
bosons in a regular fashion in order to avoid the strong re- 
pulsive interactions dominating the effective low-energy 
description. Quantum fluctuations correspond then to 
kinetic processes of the hardcore bosons in the effective 
model. Here we have treated the kinetic terms on a mean- 
field level by reducing the complex many-body problem 
to an effective one-body problem which is solved self- 
consistently. We have concentrated on crystalline solu- 
tions of the mean-field equations. 

It is the central result of this article that quantum 
fluctuations are essential to understand the low-density 
phase diagram of the compound SrCu2 (603)2- Further- 
more, we find several indications that the ratio J'/ J is 
very large J'/ J > 0.65 placing the frustrated quantum 
magnet SrCu2 (603)2 in the most complex parameter 
regime. Both aspects together result in a low-density 
phase diagram which is different to the classical sequence 



M = 1/9, M — 2/15, and M = 1/6 but much closer 
to the experimental results and their phenomenlogical 
interpretation 24 : a) We find that quantum fluctuations 
stabilize in a natural way a 1/8 plateau with a diamond 
unit cell consistent with experiments. The quantum fluc- 
tuations induced by the DM interactions play a central 
role, b) The structure of the 1/6 plateau changes al- 
ready on the classical level for J'/ J > 0.67. Micro- 
scopically, the repulsive interaction V3 becomes compa- 
rable to V4 in this parameter regime and, as a conse- 
quence, it becomes favorable to realize structures avoid- 
ing V3'. This is likely the reason why the classical struc- 
tures found in Ref. 18 for J'/ J = 0.5 arc different to the 
ones observed experimentally. All classical structures for 
J'/J = 0.5 depend strongly on V£. For J'/J > 0.67, 
one finds that the three classical structures 1/6-stripe, 
1/6-square, and 1/6- new have exactly the same classi- 
cal energy. Interestingly, the structure 1/6-stripe is pro- 
posed phenomenologically 24 . In our mean-field theory, 
quantum fluctuations favor the structure 1 / 6- new but we 
would like to stress that energy differences are very small, 
c) In the regime between 1 /8 and 1 /6 our mean-field the- 
ory predicts a sequence of stripe structures all containing 
substructures of densities 1/8 and 1/6. This sequence 
includes a 2/15 plateau. Qualitatively, this is again con- 
sistent with the phenomenological interpretation of the 
NMR data 24 . But in the mean- field treatment it is always 
the 1/8 plateau with a diamond unit cell which is realized 
inside the stripes because it benefits most from the in- 
duced quantum fluctuations, d) The mean-field solution 
is still consistent with a plateau at density 1 /9 which is in 
disagreement with the NMR data on SrCu2(BOs)2 but in 
accordance with high-field torque measurements by Se- 
bastian et al 8 . But we stress that it is well possible that 
the 1/9 plateau melts completely when superfluid solu- 
tions are considered, when effects beyond the mean-field 
level are taken into account, or when additional magnetic 
couplings are present. Indeed, as we have shown above, 
already slightly larger values for the kinetic terms lead 
to a melting of the 1 /9 plateau. 

Clearly, there are limitations of our microscopic ap- 
proach. One has to distinguish the following three levels: 
i) limitations of the mean-field treatment, ii) limitations 
of the effective low-energy description, and iii) limitations 
of the studied microscopic Hamiltonian. 

i) We have focused on crystalline solutions of the 
mean-field equations because the magnetization curve of 
SrCu2(B03)2 is dominated by magnetization plateaux 
corresponding to gapped bosonic Mott insulators. Defi- 
nitely, this approach is hard to justify in the dilute limit 
where one expects superfluid solutions to be present in 
the phase diagram. In our theory, we find a sequence of 
low-density plateau structures with a negative energy for 
jj, = fj,Q. This part of the phase diagram denoted by a in 
this work corresponds likely to a superfluid phase. Addi- 
tionally, the mean-field treatment becomes problematic 
when the physics cannot be reduced to an effective one- 
body problem. We found solutions in the density regime 
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FIG. 24: Illustrations of the relevant correlated hopping pro- 
cesses for the structures at density 1/8 (top) and density 1/6 
(down) . 



between 1/8 and 1/6 denoted by j3 where particles in 
the regions are completely delocalized between different 
dimers. This is clearly an artefact of the mean- field calcu- 
lation and we do not trust the mean-field solution in this 
density regime. One possibility to go beyond might be to 
map the effective low-energy model in terms of hardcore 
bosons to a quantum spin model which is then solved by 
spin- wave theory on various large unit cells. 

ii) The effective low-energy theory is limited by the 
extrapolations used for the various couplings and by the 
choice of operators included in the effective low-energy 
description. Qualitatively, we expect that all revelant 
tendencies are correctly described by the extrapolations 
up to the relevant regime J'/ J > 0.65. It is neverthe- 
less well possible that errors in the extrapolations lead 
to little changes, e.g. the value J'/ J > 0.67 where the 
structure of the 1/6 plateau changes already on the clas- 
sical level might shift slightly. One option to go beyond 
is to determine the effective couplings like the repulsive 
interactions non-perturbatively using continuous unitary 
transformation 40-42 or contractor renormalization 22 . In- 
deed, this would be of great help in order to describe the 
physics quantitatively in the relevant coupling regime of 
large J'/ J. The choice of operators to be treated in the 
effective low-energy models depends on the density. Here 
we have included two-body interactions and one-particle 
hopping. At low densities this is certainly a very good 
approximation. As discussed above, corrections like cor- 
related hopping are negligible for n < 1/8. We therefore 
expect that all our findings for the most relevant 1/8 
plateau are quantitatively correct. In contrast, the same 
is less obvious for the density regime 1/8 < n < 1/6. 
Here certain correlated hopping processes could be im- 
portant. A correct treatment of the correlated hopping 



processes is most likely only possible if one goes beyond a 
mean-field description. At larger densities like n = 1/4, 
correlated hopping is expected to be relevant 30 . 

iii) Our results demonstrate clearly that residual in- 
teractions of the order J/100 beyond the pure Shastry- 
Sutherland model are important for the magnetization of 
the frustrated quantum magnet SrCu2 (1303)2 • Here we 
have included the effects of DM interactions and of ad- 
ditional Heisenberg couplings to next-nearest neighbors 
by including the order-one contribution to the effective 
low-energy model. Higher orders in such couplings are 
expected to be unimportant. But let us stress that we 
have not considered effects linear in the additional cou- 
plings times a factor (J'/J) n with n > 0. Most im- 
portantly, it would be interesting to study the feedback 
effect of the additional couplings on the two-body repul- 
sive interactions 43 . Finally, let us mention that the role 
of inter-plane exchange couplings for the formation of 
plateau structures is currently not understood. 

Altogether, we have presented an important step 
towards a microscopic description of the frustrated 
quantum magnet SrCu2 (603)2- Quantum fluctuations 
induced by unfrustrated additional magnetic couplings 
like the DM interaction are very important in the low- 
density regime. A natural mechanism for the realization 
of the most prominent 1/8 plateau with a diamond 
unit cell is discovered. Furthermore, several indications 
for a rather large ratio J'/J > 0.65 are found placing 
SrCu2(B03)2 in the most challenging parameter regime. 
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VIII. APPENDIX 



A. Classical plateau energies 



In the following we list the explicit expressions for the 
classical energy per dimer e c \ = E C \/(JN) of all consid- 



ered plateaux from Fig. 4: 
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